A macroevolutionary role for chromosomal fusion and fission in Erebia butterflies.

The impact of large-scale chromosomal rearrangements, such as fusions and fissions, on speciation is a long-standing conundrum. We assessed whether bursts of change in chromosome numbers resulting from chromosomal fusion or fission are related to increased speciation rates in Erebia, one of the most species-rich and karyotypically variable butterfly groups. We established a genome-based phylogeny and used state-dependent birth-death models to infer trajectories of karyotype evolution. We demonstrated that rates of anagenetic chromosomal changes (i.e., along phylogenetic branches) exceed cladogenetic changes (i.e., at speciation events), but, when cladogenetic changes occur, they are mostly associated with chromosomal fissions rather than fusions. We found that the relative importance of fusion and fission differs among Erebia clades of different ages and that especially in younger, more karyotypically diverse clades, speciation is more frequently associated with cladogenetic chromosomal changes. Overall, our results imply that chromosomal fusions and fissions have contrasting macroevolutionary roles and that large-scale chromosomal rearrangements are associated with bursts of species diversification.

More hidden diversity in a cryptic species complex: a new subspecies of Leptideasinapis (Lepidoptera, Pieridae) from Northern Iran.

A new subspecies of Leptideasinapis from Northern Iran, discovered by means of DNA barcoding, is described as Leptideasinapistabarestanassp. nov. The new subspecies is allopatric with respect to other populations of L.sinapis and is genetically distinct, appearing as a well-supported sister clade to all other populations in COI-based phylogenetic reconstructions. Details on karyotype, genitalia, ecology and behaviour for the new subspecies are given and a biogeographical speciation scenario is proposed.

The genome sequence of the lesser marbled fritillary, Brenthis ino, and evidence for a segregating neo-Z chromosome.

The lesser marbled fritillary, Brenthis ino (Rottemburg, 1775), is a species of Palearctic butterfly. Male Brenthis ino individuals have been reported to have between 12 and 14 pairs of chromosomes, a much-reduced chromosome number than is typical in butterflies. Here, we present a chromosome-level genome assembly for Brenthis ino, as well as gene and transposable element annotations. The assembly is 411.8 Mb in length with a contig N50 of 9.6 Mb and a scaffold N50 of 29.5 Mb. We also show evidence that the male individual from which we generated HiC data was heterozygous for a neo-Z chromosome, consistent with inheriting 14 chromosomes from one parent and 13 from the other. This genome assembly will be a valuable resource for studying chromosome evolution in Lepidoptera, as well as for comparative and population genomics more generally.

Genomics Reveal Admixture and Unexpected Patterns of Diversity in a Parapatric Pair of Butterflies.

We studied the evolutionary relationship of two widely distributed parapatric butterfly species, Melitaea athalia and Melitaea celadussa, using the ddRAD sequencing approach, as well as genital morphology and mtDNA data. M. athalia was retrieved as paraphyletic with respect to M. celadussa. Several cases of mito-nuclear discordance and morpho-genetic mismatch were found in the contact zone. A strongly diverged and marginally sympatric clade of M. athalia from the Balkans was revealed. An in-depth analysis of genomic structure detected high levels of admixture between M. athalia and M. celadussa at the contact zone, though not reaching the Balkan clade. The demographic modelling of populations supported the intermediate genetic make-up of European M. athalia populations with regards to M. celadussa and the Balkan clade. However, the dissimilarity matrix of genotype data (PCoA) suggested the Balkan lineage having a genetic component that is unrelated to the athalia-celadussa group. Although narrowly sympatric, almost no signs of gene flow were found between the main M. athalia group and the Balkan clade. We propose two possible scenarios on the historical evolution of our model taxa and the role of the last glacial maximum in shaping their current distribution. Finally, we discuss the complexities regarding the taxonomic delimitation of parapatric taxa.

A DNA barcode library for the butterflies of North America.

Although the butterflies of North America have received considerable taxonomic attention, overlooked species and instances of hybridization continue to be revealed. The present study assembles a DNA barcode reference library for this fauna to identify groups whose patterns of sequence variation suggest the need for further taxonomic study. Based on 14,626 records from 814 species, DNA barcodes were obtained for 96% of the fauna. The maximum intraspecific distance averaged 1/4 the minimum distance to the nearest neighbor, producing a barcode gap in 76% of the species. Most species (80%) were monophyletic, the others were para- or polyphyletic. Although 15% of currently recognized species shared barcodes, the incidence of such taxa was far higher in regions exposed to Pleistocene glaciations than in those that were ice-free. Nearly 10% of species displayed high intraspecific variation (>2.5%), suggesting the need for further investigation to assess potential cryptic diversity. Aside from aiding the identification of all life stages of North American butterflies, the reference library has provided new perspectives on the incidence of both cryptic and potentially over-split species, setting the stage for future studies that can further explore the evolutionary dynamics of this group.

High resolution DNA barcode library for European butterflies reveals continental patterns of mitochondrial genetic diversity.

The study of global biodiversity will greatly benefit from access to comprehensive DNA barcode libraries at continental scale, but such datasets are still very rare. Here, we assemble the first high-resolution reference library for European butterflies that provides 97% taxon coverage (459 species) and 22,306 COI sequences. We estimate that we captured 62% of the total haplotype diversity and show that most species possess a few very common haplotypes and many rare ones. Specimens in the dataset have an average 95.3% probability of being correctly identified. Mitochondrial diversity displayed elevated haplotype richness in southern European refugia, establishing the generality of this key biogeographic pattern for an entire taxonomic group. Fifteen percent of the species are involved in barcode sharing, but two thirds of these cases may reflect the need for further taxonomic research. This dataset provides a unique resource for conservation and for studying evolutionary processes, cryptic species, phylogeography, and ecology.

Two ways to be endemic. Alps and Apennines are different functional refugia during climatic cycles.

Endemics co-occur because they evolved in situ and persist regionally or because they evolved ex situ and later dispersed to shared habitats, generating evolutionary or ecological endemicity centres, respectively. We investigate whether different endemicity centres can intertwine in the region ranging from Alps to Sicily, by studying their butterfly fauna. We gathered an extensive occurrence data set for butterflies of the study area (27,123 records, 269 species, in cells of 0.5 × 0.5 degrees of latitude-longitude). We applied molecular-based delimitation methods (GMYC model) to 26,557 cytochrome c oxidase subunit 1 (COI) sequences of Western Palearctic butterflies. We identified entities based on molecular delimitations and/or the checklist of European butterflies and objectively attributed occurrences to their most probable entity. We obtained a zoogeographic regionalisation based on the 69 endemics of the area. Using phylogenetic ANOVA we tested if endemics from different centres differ from each other and from nonendemics for key ecological traits and divergence time. Endemicity showed high incidence in the Alps and Southern Italy. The regionalisation separated the Alps from the Italian Peninsula and Sicily. The endemics of different centres showed a high turnover and differed in phylogenetic distances, phenology and distribution traits. Endemics are on average younger than nonendemics and the Peninsula-Sicily endemics also have lower variance in divergence than those from the Alps. The observed variation identifies Alpine endemics as paleoendemics, now occupying an ecological centre, and the Peninsula-Sicily ones as neoendemics, that diverged in the region since the Pleistocene. The results challenge the common view of the Alpine-Apennine area as a single "Italian refugium".

Incomplete Sterility of Chromosomal Hybrids: Implications for Karyotype Evolution and Homoploid Hybrid Speciation.

Heterozygotes for major chromosomal rearrangements such as fusions and fissions are expected to display a high level of sterility due to problems during meiosis. However, some species, especially plants and animals with holocentric chromosomes, are known to tolerate chromosomal heterozygosity even for multiple rearrangements. Here, we studied male meiotic chromosome behavior in four hybrid generations (F1-F4) between two chromosomal races of the Wood White butterfly Leptidea sinapis differentiated by at least 24 chromosomal fusions/fissions. Previous work showed that these hybrids were fertile, although their fertility was reduced as compared to crosses within chromosomal races. We demonstrate that (i) F1 hybrids are highly heterozygous with nearly all chromosomes participating in the formation of trivalents at the first meiotic division, and (ii) that from F1 to F4 the number of trivalents decreases and the number of bivalents increases. We argue that the observed process of chromosome sorting would, if continued, result in a new homozygous chromosomal race, i.e., in a new karyotype with intermediate chromosome number and, possibly, in a new incipient homoploid hybrid species. We also discuss the segregational model of karyotype evolution and the chromosomal model of homoploid hybrid speciation.

Rapid colour shift by reproductive character displacement in Cupido butterflies.

Reproductive character displacement occurs when competition for successful breeding imposes a divergent selection on the interacting species, causing a divergence of reproductive traits. Here, we show that a disputed butterfly taxon is actually a case of male wing colour shift, apparently produced by reproductive character displacement. Using double digest restriction-site associated DNA sequencing and mitochondrial DNA sequencing we studied four butterfly taxa of the subgenus Cupido (Lepidoptera: Lycaenidae): Cupido minimus and the taxon carswelli, both characterized by brown males and females, plus C. lorquinii and C. osiris, both with blue males and brown females. Unexpectedly, taxa carswelli and C. lorquinii were close to indistinguishable based on our genomic and mitochondrial data, despite displaying strikingly different male coloration. In addition, we report and analysed a brown male within the C. lorquinii range, which demonstrates that the brown morph occurs at very low frequency in C. lorquinii. Such evidence strongly suggests that carswelli is conspecific with C. lorquinii and represents populations with a fixed male brown colour morph. Considering that these brown populations occur in sympatry with or very close to the blue C. osiris, and that the blue C. lorquinii populations never do, we propose that the taxon carswelli could have lost the blue colour due to reproductive character displacement with C. osiris. Since male colour is important for conspecific recognition during courtship, we hypothesize that the observed colour shift may eventually trigger incipient speciation between blue and brown populations. Male colour seems to be an evolutionarily labile character in the Polyommatinae, and the mechanism described here might be at work in the wide diversification of this subfamily of butterflies.

How long is 3 km for a butterfly? Ecological constraints and functional traits explain high mitochondrial genetic diversity between Sicily and the Italian Peninsula.

Populations inhabiting Mediterranean islands often show contrasting genetic lineages, even on islands that were connected to the mainland during glacial maxima. This pattern is generated by forces acting in historical and contemporary times. Understanding these phenomena requires comparative studies integrating genetic structure, functional traits and dispersal constraints. Using as a model the butterfly species living across the Messina strait (3 km wide) separating Sicily from the Italian Peninsula, we aimed to unravel the mechanisms limiting the dispersal of matrilines and generating genetic differentiation across a narrow sea strait. We analysed the mitochondrial COI gene of 84 butterfly species out of 90 documented in Sicily and compared them with populations from the neighbouring southern Italian Peninsula (1,398 sequences) and from the entire Palaearctic region (8,093 sequences). For each species, we regressed 13 functional traits and 2 ecological constraints to dispersal (winds experienced at the strait and climatic suitability) against genetic differentiation between Sicily and Italian Peninsula to understand the factors limiting dispersal. More than a third of the species showed different haplogroups across the strait and most of them also represented endemic haplogroups for this island. One fifth of Sicilian populations (and 32.3% of endemic lineages) had their closest relatives in distant areas, instead of the neighbouring Italian Peninsula, which suggests high relictuality. Haplotype diversity was significantly explained by the length of the flight period, an intrinsic phenology trait, while genetic differentiation was explained by both intrinsic traits (wingspan and degree of generalism) and contemporary local constraints (winds experienced at the strait and climatic suitability). A relatively narrow sea strait can produce considerable differentiation among butterfly matrilines and this phenomenon showed a largely deterministic fingerprint. Because of unfavourable winds, populations of the less dispersive Sicilian butterflies tended to differentiate into endemic variants or to maintain relict populations. Understanding these phenomena required the integration of DNA sequences, species traits and physical constraints for a large taxon at continental scale. Future studies may reveal if the patterns here shown for mitochondrial DNA are also reflected in the nuclear genome or, alternatively, are the product of limited female dispersal.

Dissecting the Effects of Selection and Mutation on Genetic Diversity in Three Wood White (Leptidea) Butterfly Species.

The relative role of natural selection and genetic drift in evolution is a major topic of debate in evolutionary biology. Most knowledge spring from a small group of organisms and originate from before it was possible to generate genome-wide data on genetic variation. Hence, it is necessary to extend to a larger number of taxonomic groups, descriptive and hypothesis-based research aiming at understanding the proximate and ultimate mechanisms underlying both levels of genetic polymorphism and the efficiency of natural selection. In this study, we used data from 60 whole-genome resequenced individuals of three cryptic butterfly species (Leptidea sp.), together with novel gene annotation information and population recombination data. We characterized the overall prevalence of natural selection and investigated the effects of mutation and linked selection on regional variation in nucleotide diversity. Our analyses showed that genome-wide diversity and rate of adaptive substitutions were comparatively low, whereas nonsynonymous to synonymous polymorphism and substitution levels were comparatively high in Leptidea, suggesting small long-term effective population sizes. Still, negative selection on linked sites (background selection) has resulted in reduced nucleotide diversity in regions with relatively high gene density and low recombination rate. We also found a significant effect of mutation rate variation on levels of polymorphism. Finally, there were considerable population differences in levels of genetic diversity and pervasiveness of selection against slightly deleterious alleles, in line with expectations from differences in estimated effective population sizes.

The conundrum of species delimitation: a genomic perspective on a mitogenetically super-variable butterfly.

The Palaearctic butterfly Melitaea didyma stands out as one of the most striking cases of intraspecific genetic differentiation detected in Lepidoptera: 11 partially sympatric mitochondrial lineages have been reported, displaying levels of divergence of up to 7.4%. To better understand the evolutionary processes underlying the diversity observed in mtDNA, we compared mtDNA and genome-wide SNP data using double-digest restriction site-associated DNA sequencing (ddRADseq) results from 93 specimens of M. didyma ranging from Morocco to eastern Kazakhstan. We found that, between ddRADseq and mtDNA results, there is a match only in populations that probably remained allopatric for long periods of time. Other mtDNA lineages may have resulted from introgression events and were probably affected by Wolbachia infection. The five main ddRADseq clades supported by STRUCTURE were parapatric or allopatric and showed high pairwise FST values, but some were also estimated to display various levels of gene flow. Melitaea didyma represents one of the first cases of deep mtDNA splits among European butterflies assessed by a genome-wide DNA analysis and reveals that the interpretation of patterns remains challenging even when a high amount of genomic data is available. These findings actualize the ongoing debate of species delimitation in allopatry, an issue probably of relevance to a significant proportion of global biodiversity.

Global invasion history of the agricultural pest butterfly Pieris rapae revealed with genomics and citizen science.

The small cabbage white butterfly, Pieris rapae, is a major agricultural pest of cruciferous crops and has been introduced to every continent except South America and Antarctica as a result of human activities. In an effort to reconstruct the near-global invasion history of P. rapae, we developed a citizen science project, the "Pieris Project," and successfully amassed thousands of specimens from 32 countries worldwide. We then generated and analyzed nuclear (double-digest restriction site-associated DNA fragment procedure [ddRAD]) and mitochondrial DNA sequence data for these samples to reconstruct and compare different global invasion history scenarios. Our results bolster historical accounts of the global spread and timing of P. rapae introductions. We provide molecular evidence supporting the hypothesis that the ongoing divergence of the European and Asian subspecies of P. rapae (∼1,200 y B.P.) coincides with the diversification of brassicaceous crops and the development of human trade routes such as the Silk Route (Silk Road). The further spread of P. rapae over the last ∼160 y was facilitated by human movement and trade, resulting in an almost linear series of at least 4 founding events, with each introduced population going through a severe bottleneck and serving as the source for the next introduction. Management efforts of this agricultural pest may need to consider the current existence of multiple genetically distinct populations. Finally, the international success of the Pieris Project demonstrates the power of the public to aid scientists in collections-based research addressing important questions in invasion biology, and in ecology and evolutionary biology more broadly.

Lack of gene flow: Narrow and dispersed differentiation islands in a triplet of Leptidea butterfly species.

Genome scans in recently separated species can inform on molecular mechanisms and evolutionary processes driving divergence. Large-scale polymorphism data from multiple species pairs are also key to investigate the repeatability of divergence-whether radiations tend to show parallel responses to similar selection pressures and/or underlying molecular forces. Here, we used whole-genome resequencing data from six wood white (Leptidea sp.) butterfly populations, representing three closely related species with karyomorph variation, to infer the species' demographic history and characterize patterns of genomic diversity and differentiation. The analyses supported previously established species relationships, and there was no evidence for postdivergence gene flow. We identified significant intraspecific genetic structure, in particular between karyomorph extremes in the wood white (L. sinapis)-a species with a remarkable chromosome number cline across the distribution range. The genomic landscapes of differentiation were erratic, and outlier regions were narrow and dispersed. Highly differentiated (FST ) regions generally had low genetic diversity (θπ ), but increased absolute divergence (DXY ) and excess of rare frequency variants (low Tajima's D). A minority of differentiation peaks were shared across species and population comparisons. However, highly differentiated regions contained genes with overrepresented functions related to metabolism, response to stimulus and cellular processes, indicating recurrent directional selection on a specific set of traits in all comparisons. In contrast to the majority of genome scans in recently diverged lineages, our data suggest that divergence landscapes in Leptidea have been shaped by directional selection and genetic drift rather than stable recombination landscapes and/or introgression.

Integrating three comprehensive data sets shows that mitochondrial DNA variation is linked to species traits and paleogeographic events in European butterflies.

Understanding the dynamics of biodiversity, including the spatial distribution of genetic diversity, is critical for predicting responses to environmental changes, as well as for effective conservation measures. This task requires tracking changes in biodiversity at large spatial scales and correlating with species functional traits. We provide three comprehensive resources to understand the determinants for mitochondrial DNA differentiation represented by (a) 15,609 COI sequences and (b) 14 traits belonging to 307 butterfly species occurring in Western-Central Europe and (c) the first multi-locus phylogenetic tree of all European butterfly species. By applying phylogenetic regressions we show that mitochondrial DNA spatial differentiation (as measured with GST , G'ST , D and DST ) is negatively correlated with species traits determining dispersal capability and colonization ability. Thanks to the high spatial resolution of the COI data, we also provide the first zoogeographic regionalization maps based on intraspecific genetic variation. The overall pattern obtained by averaging the spatial differentiation of all Western-Central European butterflies shows that the paradigm of long-term glacial isolation followed by rapid pulses of post-glacial expansion has been a pervasive phenomenon in European butterflies. The results and the extensive data sets we provide here constitute the basis for genetically-informed conservation plans for a charismatic group in a continent where flying insects are under alarming decline.

A mirage of cryptic species: Genomics uncover striking mitonuclear discordance in the butterfly Thymelicus sylvestris.

Mitochondrial DNA (mtDNA) sequencing has led to an unprecedented rise in the identification of cryptic species. However, it is widely acknowledged that nuclear DNA (nuDNA) sequence data are also necessary to properly define species boundaries. Next generation sequencing techniques provide a wealth of nuclear genomic data, which can be used to ascertain both the evolutionary history and taxonomic status of putative cryptic species. Here, we focus on the intriguing case of the butterfly Thymelicus sylvestris (Lepidoptera: Hesperiidae). We identified six deeply diverged mitochondrial lineages; three distributed all across Europe and found in sympatry, suggesting a potential case of cryptic species. We then sequenced these six lineages using double-digest restriction-site associated DNA sequencing (ddRADseq). Nuclear genomic loci contradicted mtDNA patterns and genotypes generally clustered according to geography, i.e., a pattern expected under the assumption of postglacial recolonization from different refugia. Further analyses indicated that this strong mtDNA/nuDNA discrepancy cannot be explained by incomplete lineage sorting, sex-biased asymmetries, NUMTs, natural selection, introgression or Wolbachia-mediated genetic sweeps. We suggest that this mitonuclear discordance was caused by long periods of geographic isolation followed by range expansions, homogenizing the nuclear but not the mitochondrial genome. These results highlight T. sylvestris as a potential case of multiple despeciation and/or lineage fusion events. We finally argue, since mtDNA and nuDNA do not necessarily follow the same mechanisms of evolution, their respective evolutionary history reflects complementary aspects of past demographic and biogeographic events.

Towards a Sedimentology of Information Infrastructures: a Geological Approach for Understanding the City.

Drawing primarily on ethnographic research performed in a city in Romania, this paper provides a thick description of police practices and information systems in that municipality. It shows various ways in which technologies mediate policing practitioners' perceptions, decisions and actions. Bringing some additional material from a case in the Dutch police in which they build risk profiles predicated on real-time data from a sensor network, the paper highlights new phenomena with ethical implications emerging at the intersection of information infrastructures and policing practices. In particular, it shows a solidifying effect of technologically mediated suspicion, a sedimentation process in technological infrastructures and the formation of pockets of prejudice in the layers of software code. To adequately account for these phenomena and others, the paper provides a set of arguments and conditions for a sedimentology of infrastructures. At the same time, it offers the first steps in a larger research project that would adapt and test the limits of a geological vocabulary and approach to understand smart urban environments.

Versatility of multivalent orientation, inverted meiosis, and rescued fitness in holocentric chromosomal hybrids.

Chromosomal rearrangements (e.g., fusions/fissions) have the potential to drive speciation. However, their accumulation in a population is generally viewed as unlikely, because chromosomal heterozygosity should lead to meiotic problems and aneuploid gametes. Canonical meiosis involves segregation of homologous chromosomes in meiosis I and sister chromatid segregation during meiosis II. In organisms with holocentric chromosomes, which are characterized by kinetic activity distributed along almost the entire chromosome length, this order may be inverted depending on their metaphase I orientation. Here we analyzed the evolutionary role of this intrinsic versatility of holocentric chromosomes, which is not available to monocentric ones, by studying F1 to F4 hybrids between two chromosomal races of the Wood White butterfly (Leptidea sinapis), separated by at least 24 chromosomal fusions/fissions. We found that these chromosomal rearrangements resulted in multiple meiotic multivalents, and, contrary to the theoretical prediction, the hybrids displayed relatively high reproductive fitness (42% of that of the control lines) and regular behavior of meiotic chromosomes. In the hybrids, we also discovered inverted meiosis, in which the first and critical stage of chromosome number reduction was replaced by the less risky stage of sister chromatid separation. We hypothesize that the ability to invert the order of the main meiotic events facilitates proper chromosome segregation and hence rescues fertility and viability in chromosomal hybrids, potentially promoting dynamic karyotype evolution and chromosomal speciation.

Use of genetic, climatic, and microbiological data to inform reintroduction of a regionally extinct butterfly.

Species reintroductions are increasingly used as means of mitigating biodiversity loss. Besides habitat quality at the site targeted for reintroduction, the choice of source population can be critical for success. The butterfly Melanargia russiae (Esper´s marbled white) was extirpated from Hungary over 100 years ago, and a reintroduction program has recently been approved. We used museum specimens of this butterfly, mitochondrial DNA data (mtDNA), endosymbiont screening, and climatic-similarity analyses to determine which extant populations should be used for its reintroduction. The species displayed 2 main mtDNA lineages across its range: 1 restricted to Iberia and southern France (Iberian lineage) and another found throughout the rest of its range (Eurasian lineage). These 2 lineages possessed highly divergent wsp alleles of the bacterial endosymbiont Wolbachia. The century-old Hungarian specimens represented an endemic haplotype belonging to the Eurasian lineage, differing by one mutation from the Balkan and eastern European populations. The Hungarian populations of M. russiae occurred in areas with a colder and drier climate relative to most sites with extant known populations. Our results suggest the populations used for reintroduction to Hungary should belong to the Eurasian lineage, preferably from eastern Ukraine (genetically close and living in areas with the highest climatic similarity). Materials stored in museum collections can provide unique opportunities to document historical genetic diversity and help direct conservation.

Gene expression profiling across ontogenetic stages in the wood white (Leptidea sinapis) reveals pathways linked to butterfly diapause regulation.

In temperate latitudes, many insects enter diapause (dormancy) during the cold season, a period during which developmental processes come to a standstill. The wood white (Leptidea sinapis) is a butterfly species distributed across western Eurasia that shows photoperiod-induced diapause with variation in critical day-length across populations at different latitudes. We assembled transcriptomes and estimated gene expression levels at different developmental stages in experimentally induced directly developing and diapausing cohorts of a single Swedish population of L. sinapis to investigate the regulatory mechanisms underpinning diapause initiation. Different day lengths resulted in expression changes of developmental genes and affected the rate of accumulation of signal molecules, suggesting that diapause induction might be controlled by increased activity of monoamine neurotransmitters in larvae reared under short-day light conditions. Expression differences between light treatment groups of two monoamine regulator genes (DDC and ST) were observed already in instar III larvae. Once developmental pathways were irreversibly set at instar V, a handful of genes related to dopamine production were differentially expressed leading to a significant decrease in expression of global metabolic genes and increase in expression of genes related to fatty acid synthesis and sequestration. This is in line with a time-dependent (hour-glass) model of diapause regulation where a gradual shift in the concentration of monoamine neurotransmitters and their metabolites during development of larvae under short-day conditions leads to increased storage of fat, decreased energy expenditures, and ultimately developmental stasis at the pupal stage.